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^ , Abstract. We show that well-chosen Lagrangians for a class of two-dimensional 

' integrable lattice equations obey a closure relation when embedded in a higher 

dimensional lattice. On the basis of this property we formulate a Lagrangian 
| description for such systems in terms of Lagrangian multiforms. We discuss the 

connection of this formalism with the notion of multidimensional consistency, and 
the role of the lattice from the point of view of the relevant variational principle. 

d 

1. Introduction 

In recent years there has been a growing interest in the integrability of discrete systems 
defined on two- or multidimensional lattices. In part the study of such systems may be 

00 . 

motivated by the search for accurate approximations to continuous systems. However, 
the modern point of view is that such lattice systems are important in their own right 
from a theoretical perspective, and, in fact, are thought to be richer and more generic 
Q\ ', than their continuous counterparts. Discrete systems have also been proposed in physics 
to model the fundamental interactions on the scale of the Planck constant where space 



and time themselves can be thought of as being discrete [201 139] - 

The earliest examples of integrable lattice systems go back to the mid 1970s and 
early 1980s, when the research was focused on discretizing known continuous soliton 
systems [U El [HI HSJ [TOl EH EH]- In recent years, the insight has developed that 
the key aspect of integrability resides in the property of multidimensional consistency 
[321 El [33]. This property entails that an equation can be embedded in a consistent way 
in a multidimensional lattice, i.e. by imposing a copy of the equation (with appropriate 
lattice parameters) in each pair of directions, such that there is no inconsistency or 
multivaluedness occurring in the evaluation of the dependent variables on each lattice 
site. Using this property a classification of two-dimensional scalar integrable lattice 
systems was given, in the affine linear case, by Adler, Bobenko and Suris 0, H] (resulting 
in what is hereafter referred to as the ABS list). In addition to the known examples 
of lattice systems of KdV type, cf. [HI [151 EH EH], this provides us with some new 
examples of integrable scalar lattice equations. 

Conventionally the Hamiltonian has been the central object in (continuous) 
integrable systems [12]. Of course, it is often possible to pass between Lagrangian 
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and Hamiltonian theories via Legendre transforms, although this is in many (non- 
Newtonian) cases not a trivial matter. Nevertheless, most integrable partial differential 
equations seem to admit a Lagrangian description; in fact, a universal Lagrangian 
structure for integrable systems admitting a Lax pair was formulated by Zakharov 
and Mikhailov |43j. In the discrete case one can argue that the Lagrangian is the 
more fundamental object^, and in fact Lagrangian structures have been established for 
several discrete integrable systems such as Lagrangian mappings [HI l4*2l [26] . These 
Lagrangian descriptions are based on a discrete calculus of variations as developed 
earlier by CadzowJT], Logan[23] and Maeda[2l] outside the scope of integrable systems. 
Furthermore, Lagrangians and/or actions were also constructed for integrable two- 
dimensional lattice equations, cf. [11 [291 [3]. Discrete Lagrangian systems on arbitrary 
graphs were proposed in [34], and a discrete variational complex was set up in [16J. 

The usual point of view is that the Lagrangian is a scalar object (or equivalently 
a volume form), which through the Euler-Lagrange equations provides us with one 
single equation (i.e. one per component of the dependent variable). In contrast, we 
take the point of view that in the case of an integrable system, where due to the 
multidimensional consistency several equations can be imposed simultaneously on one 
and the same dependent variable, the Lagrangian should actually be an extended object 
capable of producing a multitude of consistent equations from a variational principle. 
Thus we propose in this paper an action in which the key ingredient is a Lagrangian 
2-form (in the case of integrable discrete equations in two independent variables) or, 
more generally, a multiform (in the case of a larger number of independent variables). 
Although the notion of a Lagrangian multiform is not new, and goes back to Cartan and 
Lepage [HI [21], cf. also [17] for a review, even in those theories the role of the Lagrangian 
is that of a volume form producing the equations of motion in a conventional way. 

In the original ABS paper [3], action functionals were given for the whole list 
of lattice equations in their classification. In the present paper we reformulate the 
Lagrangian structures of these lattice systems, identifying a specific form of the 
Lagrangians, and we show by explicit computation that for each case considered, a 
closure-type relation holds when we embed these systems into a higher dimensional 
lattice. This closure relation effectively indicates that the discrete Lagrangian is a 
closed 2-form on the multidimensional lattice, and we consider this as a manifestation 
of multidimensional consistency on the Lagrangian level. We argue that since this 
relation implies surface independence of the relevant action, the variational principle 
that describes multidimensionally consistent systems requires variations with respect 
not only to the field variables, but also to the geometry of the lattice on which the 
integrable discrete equation is defined. By presenting an explicit continuous example 
of a multidimensionally consistent system, namely the generating PDE associated with 
the lattice KdV equation [30] , we draw in the penultimate section some analogies with 
the continuous case. 

% Drawing the parallel with quantum mechanics, we adhere to Dirac's opinion [llj . 
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2. Multidimensional Consistency of Lattices and 3-point Lagrangians 

In this paper, following [3], we consider lattice systems of the following form: there 
are two independent (discrete) variables ni,n 2 corresponding to two lattice directions, 
two lattice parameters a%, a 2 (which can be thought of as measures for the grid size) 
associated with the n±, n 2 directions respectively, and a scalar dependent variable 
u{n\,ri2)- The parameters a%, a 2 may take on continuous values, whilst n\,n 2 are 
discrete coordinates of the lattice, but may in principle take on continuous values as 
long as elementary shifts in these variables amount to increments by one unit. For 
ease of notation, let u = u(ni,n 2 ), u\ = u(n\ + l,n 2 ) and u 2 = u(ni,n 2 + 1), as 
indicated in Figured! Backwards shifts in u are denoted by = u{rii — l,n 2 ) and 
w_2 = u(ni,n 2 — 1). In [3] the classification problem of quadrilateral lattice equations 
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Figure 1. 2-d lattice 

of the following form was considered 

Q(u,ux, u 2 ,u h2 ; a>i,a 2 ) = (2.1) 

where Q is affine linear and subject to the symmetry of the square (D4 symmetry). The 
main classification criterion is that of multidimensional consistency, by which we mean 
the following. 

It is known [32} [5] that these equations can be consistently embedded in a 
3-dimensional lattice, imposing the same form of the equation (with appropriate 
parameters) in all 2-dimensional sublattices. This has appeared in many places 
in the recent literature, but for the sake of self-containedness we present it again 
here. Imposing in addition to the equation Q(u, Ui, u 2 , 1*1,2; «i, ol 2 ) = the equations 
Q(u, u 2 , U3, 1*2,3; «2, CH3) = and Q(u, U3, u\, 1*1,3; 0^3, ol\) = on the other faces of the 
elementary cube of Figure [21 and given the initial values u, u±, u 2 , m 3 , there are in 
principle three ways in which to compute the value of Wi,2,3- If these three values 
coincide, the equation is called consistent-around-the-cube, and is considered to be 
multi-dimensionally consistent. As was pointed out in [3Tj this phenomenon is analogous 
to the existence of integrable hierarchies in nonlinear evolution equations of soliton type. 
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Figure 2. Elementary cube 



The result of the classification study of [3] was a list of 9 equations, up to Mobius 
transformations, labelled H1-H3, Q1-Q4 and A1-A2. Some of these were already well- 
known, e.g. the discrete KdV equation [Hj, but the classification also produced some 
new equations. Furthermore, it was shown in [3] that all these equations admit an action 
principle, which is based on the so-called 3-leg form of the quadrilateral equation. In 
some cases, namely lattice equations "of KdV type" a Lagrangian description had been 
previously established, (SJ 29j. From the actions given in [3], one can infer 4-point 
Lagrangians, however, for our purpose it is more useful to identify 3-point Lagrangians. 
In terms of these the action will take the form 



and in this specific form Lagrangians of all ABS equations can be established. 

The discrete Euler-Lagrange equations arising from the variational principle that 
5S = 0, under local variations 5u{ni, n 2 ) of the dependent variable, are given by 



Cj(u, ui, u 2 ; ai, a 2 ) + £>(u-i,u, u-i j2 ] «i, a 2 ) + £(u_ 2 , mi,-2> w, «i, oc 2 ) =0 (2.3) 



Below we list specific examples of ABS lattice equations together with their 3-point 
Lagrangians. Although similar formulae can be established for the remaining cases 
in the ABS list, we will restrict ourselves here to these particular examples for the 
remainder of the paper. It should be noted that the discrete Euler-Lagrange equations 
( 12. 3j) do not give the quadrilateral lattice equations themselves, but rather a discrete 
derivative of the original equation which is defined on 7 points of the lattice (lattice 
equations on 7-point stencils have attracted a considerable amount of interest in recent 
years, cf. e.g. [27J). The Euler-Lagrange equation actually results in a compound of 
two copies of the 3-leg form of the original equation: one reflected in the r^-direction 
and one reflected in the ^-directional 

§ We note that in the generic case, the D4 symmetry is not always manifest on the level of the 3-leg 




(2.2) 
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2.1. HI 

This is the discrete potential Korteweg de Vries equation, one of the most fundamental 
examples in discrete integrable systems. The Lagrangian was first given in [8]. 
The original equation is 

(u — iti j2 ) (iti — u 2 ) — ai + a 2 — 0; (2.4a) 
written in 3-leg form this is 

(if + ui) - (u + u 2 ) = ai ~ a2 (2.4b) 

U - U 1>2 

and it possesses the Lagrangian 

£ = {u\ — u 2 )u — («x — a 2 ) ln(-ux — u 2 ) (2.4c) 

which through the Euler-Lagrange equation (12.31) leads to the 7-point equation 

di — a 2 ot\ — a 2 
Mi — m_ 2 1- W-i ~ u 2 = (2.4d) 

which contains two copies of the 3-leg form ( I2.4bl) . We get a similar 7-point equation, 
consisting of two copies of the original 4-point equation, from all of the Lagrangians 
given below. 

2.2. H2 

The original equation is 

(u - Ui )2 )(iii - u 2 ) - («i - a 2 )(u + ui + u 2 + u 1>2 ) - a\ + a\ = 0; (2.5a) 

written in 3-leg form this is 

u + u x + Oi _ m - Mi ]2 + «i - a 2 ^ 

"U + U 2 + « 2 M — Ml, 2 — Oil + «2 

and it possesses the Lagrangian 

£ = (u + Mi + «i) ln(u + Mi + cti) - (m + m 2 + a 2 ) ln(u + u 2 + a 2 ) 
— (mi — m 2 + «i — a 2 ) ln(Mi — m 2 + a x — a 2 ) 

+ (ui — u 2 — ax + Oi 2 ) \n(ui — u 2 — ai + Oi 2 ). (2.5c) 

form, e.g. in the case of Q4 where the 3-leg form lives on the level of the uniformizing variables of the 
relevant elliptic curve. It is through the connection with the affine linear form of the equations that 
the symmetry under reversal of the shifts becomes apparent. 
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This is also known as the discrete modified (potential) KdV equation. 
The original equation is 

Ui(uui + u 2 Ui j2 ) — ct 2 (uu 2 + UxUip) + 8(a\ — a 2 ) = 0; 

written in 3-leg form this is 

e 2x+2x! _|_ ^ g 2ai s inh(x — X\ j2 — CL\ + a 2 ) 



e 2x+2x 2 _|_ ^ e 2a2 sinh(x — xi j2 + ai — a 2 ) 

where it = e 2x and «i = e 2ai , and it possesses the Lagrangian 

. T . , uui . uu 2 ,a 2 ui. ,aiux, 

£ = - Li 2 ( =) + Li 2 ( -) + Li 2 ( ) - Li 2 ( 

— aid —a 2 o a\u 2 a 2 u 2 



+ ln(«)+ln(a 2 2 )ln^ 



where Li 2 (z) is the dilogarithm function given in flA.ip of appendix A, where 
identities for the dilogarithm function are presented. 

2.1 Ql\s=o 

The original equation is 

ai(u - u 2 )(ui - ui )2 ) - ol 2 {u - ui)(u 2 - it 1>2 ) = 0; 
written in 3-leg form this is 

«i a 2 ct\ — a 2 

U — Ui U — U 2 U — U\ 2 

and it possesses the Lagrangian 

L = a.\ ln(it — U\) — a 2 ln(u — u 2 ) — (pt\ — a 2 ) ln(iti — u 2 ). 

2.5. Ql\ m 

The original equation is 

ai(u - u 2 )(ui - «i >2 ) - a 2 (u - Ui)(u 2 - u 1;2 ) + 5 2 a 1 a 2 (a 1 - a 2 ) = 0; 
written in 3-leg form this is 

u — u\ + a>i5\ I u — u 2 — a 2 5\ f u ~ u i,2 + oti8 — a 2 S 



v m — u\ — a\5 J \u — u 2 + a 2 5 J \u — iti j2 — a,\8 + a 2 5 
and it possesses the Lagrangian 

£ = (« — «! + a\8) In (it — Mi + a±8) — (it — u\ — a±8) ln(it — iti — a\8) 

— (u — it 2 + a 2 <$) ln(ix — it 2 + a 2 S) + (it — u 2 — a 2 5) ln(it — it 2 — a 2 5) 

— (u\ — u 2 + q;i5 — a 2 5) ln(iti — it 2 + ai<5 — a 2 5) 
+ («! — u 2 — a.\8 + a 2 ^) ln(iii — u 2 — a.\8 + a 2 5). 
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Written in a slightly different form, this equation is known as the Homotopy equation, 
and appears in the literature in [28]. The original equation is 

(a\ - a\) (mti )2 + UxU 2 ) + a 2 (a\ - 1) (uux + u 2 u lt2 ) - ctx{oc\ - 1) (uu 2 + MiWi, 2 ) = 0; (2.9a) 

written in 3-leg form this is 

'sinh(x — xi + ai)\ / sinh(a; — x 2 — a 2 )\ f sinh(a; — x 1)2 + ax — a 2 ) 



v sinh(x — x\ — a\) J \sinh(x — x 2 + a 2 ) J \sinh(x — x\ >2 , — a\ + a 2 

(2.9b) 

where u = e 2x and a,\ = e 2ai , and it possesses the Lagrangian 

z=-u 2 [^)+u 2 (^) + uJ^)-uJ u 



u\ j \ol\U\) \ u 2 j \a 2 u 2 
+ Li 2 ( ™ )- Li 2 ( )+ l n(a? ) J2*±). (2.9c) 

2.7. Al 

The original equation is 

a x {u + u 2 )(ui + u lj2 ) - a 2 (u + Ui)(u 2 + u 1:2 ) - 5 2 a 1 a 2 (a 1 - a 2 ) = 0; (2.10a) 

written in 3-leg form this is 

/ u + Mi + a\5\ ( u + u 2 - a 2 5 \ _ ( u - it 1>2 + a x 8 - a 2 S \ 
\u + u\ — ai5 J \u + u 2 + a 2 5 J \u — «i >2 — ot\8 + a 2 5 J 

and it possesses the Lagrangian 

L = {u + ui + ot\8) In (it + U\ + cx\5) — (u + Ux — aiS) ln(it + iti — «i5) 

- (it + u 2 + a 2 8) ln(it + u 2 + a 2 5) + (it + u 2 - a 2 5) ln(it + w 2 — a 2 5) 

— (w 2 — iii + ol\5 — a 2 5) ln(u 2 — iii + ol\5 — a 2 <5) 

+ (it 2 — iti — ot\8 + a 2 S) ln(it 2 — iti — cti# + a 2 5). (2.10c) 

The original equation is 

(al — al)(uuiu 2 ui t2 + l) + a 2 (al — l)(uu 2 + uiUi t2 )—ai(al — l)(uui + u 2 ui j2 ) = 0; (2.11a) 

written in 3-leg form this is 

/ sinh(x + X\ + a{) \ ( sinh(x + x 2 — a 2 ) \ / sinh(x — xx,2 + a\ — a 2 ) \ 
\ v sinh(x + x\ — ax) J \sinh(x + x 2 + a 2 ) J \sxnh.{x — xx,2 — «i + a 2 ) J 

(2.11b) 
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where u = e 2x and ct\ = e 2ai , and it possesses the Lagrangian 

£ = - Li 2 («i-uwi) + Li 2 ( — ^ ) + Li 2 (a 2 uu 2 ) - Li 2 

V a x J 





) 



+ Li 2 



a 2 ui 



OL\U 2 



) 



- Li 2 



) 



+ ln(aj) In 



a 2 u 2 



(2.11c) 



For the Lagrangians of the cases given in the list 2.1-2.8, we will next establish an 
important new property. 

3. Closure relation and Lagrangian 2-forms 

The main observation of this paper is that all the lattice systems, together with their 
3-point Lagrangians, as given in the previous section, possess a remarkable property 
which we refer to as the closure relation, when we embed both the equation and the 
Lagrangian in a 3-dimensional lattice. In order to formulate this property we introduce 
the notation of the difference operator Aj which acts on functions / of u = u(ni, n 2 , n 3 ) 
by the formula Aif(u) = f(ui) — f{u), and on a function g of u and its shifts by the 
formula Aig(u,Uj,Uk) = g(ui,u it j,u ijk ) — g(u,Uj,u k ), in which, as before, the suffix i 
denotes a shift in the direction associated with the variable n^. The following statement 
holds true. 

Proposition: 

All the 3-point Lagrangians given in the list 2.1-2.8 when embedded in a three- 
dimensional lattice, as explained in section 2, satisfy the following relation on solutions 
of the quadrilateral lattice system: 

AiL(u, it 2 , u 3 ; a 2 , a 3 ) + A 2 L(u, u 3 , ux] a 3 , a x ) + A 3 £j(u, m, u 2 ; a 1} a 2 ) = (3.1) 

This can be established by explicit computation, and has been verified in all cases 
in the list 2.1-2.8. Below we will demonstrate this computation for the case of HI. 
Furthermore, in appendix B we will present the computation in the case of H3, which is 
somewhat more involved and relies on a number of identities for the dilogarithm function 
Li 2 , see e.g. [221 EE], the relevant ones of which have been reproduced in appendix A. 
For the Lagrangians of the remaining equations in the ABS list the computations are 
more implicit and we delegate those to a future publication. 

Example: HI To illustrate the proposition in the simplest case, we perform the 
following computation. By definition of the Lagrangians we have 

Ai£(tt, u 2 , u 3 ; a 2 , a 3 ) + A 2 L(u, u 3 , u x ; a 3 , a x ) + A 3 L(u, u x ,u 2 \ a lf a 2 ) 

= (wi,2 - u^ux - (a 2 - a 3 ) ln(u lt2 - ux, 3 ) - [u 2 - u 3 )u 

+ («2 - a 3 ) \rx(u 2 - u 3 ) + (u 2 ,3 - "1,2)^2 - (a 3 - ax) ln(u 2 , 3 - u lt2 ) 

- (u 3 - ux)u + (a 3 - ax) \n(u 3 - ux) + (u lj3 - u 2t3 )u 3 

- (ax - a 2 ) \n(u 1:3 - u 2j3 ) - (ux - u 2 )u + (ax - a 2 ) \n(ux - u 2 ). (3.2) 
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(3.3) 



where A lj2)3 is invariant under permutations of the indices, the expression (13. 2p reduces 
to 



+M,2, 3 {u 3 - Ui)u 2 - ("3 - aii) hi(A 1)2 , 3 (M3 - «i)) - («3 - + (a 3 - ai) ln(« 3 - «i) 



where we have tried to organize the succession of terms to make it manifest which 
groupings of terms cancel out against each other. 

In a similar way the cases of H2, Ql and Al can be verified, whereas the cases of 
H3, Q3|<5 = o and A2 all involve the dilogarithm function Li 2 and can be verified along 
similar lines as the computation in appendix B. 

In order to discuss the implications of the statement above, we need to introduce 
some further notation. Let denote the unit vector in the lattice direction labelled by 
i and let any point in the multidimensional lattice be specified by the vector n whose 
components are the coordinates n±, n 2 , . . . of the latticqjj], then elementary shifts in the 
lattice can be generated by the action n — > n + e^. Specifying an elementary oriented 
plaquette in this lattice requires the following data: the position n of one of its vertices 
in the lattice and the lattice directions given by the base vectors e^e^. One way to 
characterize the oriented plaquette is by the ordered triplet o~ij(n) = (n, n + ei, n + ej). 
Since the 3-point Lagrangians depend on two directions in the lattice, and when 
embedded in a multidimensional lattice at each point can be associated with an oriented 
plaquette ^(n), we can think of these Lagrangians as defining a discrete 2- form ^(n) 
whose evaluation on that plaquette is given by the Lagrangian function as follows 



Here it is understood that (I3.5P is the contribution to the action functional associated 
with the given plaquette o~ij{n) as described above. Choosing now a surface a in the 
multidimensional lattice consisting of a connected configuration of elementary plaquettes 
o~ij(n), such as illustrated in Figured (which could be an infinite surface or a compact 
surface, with or without boundary) we can define an action on that surface simply by 

|| If we select a lattice of finite dimensionality we could write the coordinates on the lattice as 
n = (m, rig, . . . ), where the lattice directions are labelled according to the natural numbers. However in 
principle one could also have an infinite dimensional lattice and even a lattice labelled by an uncountable 
set. 



Ai, 2j3 (u 2 - u 3 )u! - (a 2 - a 3 ) ln^4i )2 ,3(«2 - u 3 )) - (u 2 - u 3 )u + (a 2 - a 3 ) ln(w 2 - u 3 ) 




Lij{n) = £j(u(n),u(n + ei),u(n + e,-); a^, ctj). 



(3.5) 
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Figure 3. Example of a surface with boundary 

summing up the contributions from each of the plaquettes on the surface, taking 
into account the directions associated with each face, i.e. we perform the sum: 

S = S[u{n);a] = J2^ = M n )- ( 3 - 6 ) 

The sum in (13.61) is unambiguous for two reasons: first, because all the Lagrangians 
considered in 2.1-2.8 have the property of antisymmetry up to a constant with respect 
to transformations i <-> j, i.e. £jij(n) = —Lji(n) + constant; second, we choose the base 
point n in such a way that <Cy(n), defined on o~ij(n), involves u(n) along with its shifts 
only in the positive i and j directions. We choose throughout this paper not to use the 
abstract notation of difference forms, cf. e.g. [25], because we want to demonstrate on 
the basis of the examples given that all statements can be established through concrete 
computations. 

It is obvious from (13.61) that the geometry of the surface a forms an integral part 
of the action functional. The closure relation (13.11) implies the invariance of the action 
under local deformations of the surface a while fixing its boundary. This we can easily 
see by considering an elementary variation of a locally flat surface at a single plaquette 
illustrated by Figure HI If 5" is the value of the action functional for the undeformed 




Figure 4. Local deformation of a discrete surface a to a surface a' 

surface in Figure H] the value for the deformed surface in Figure H] can be computed as 
follows 

S' = S - Cj(u, u h Uj] aci, oij) + £(itfc, u i>k , u j)k ; a it a,-) + u id , u i>k ; a>j, a k ) 

+ L(uj,u jjk , u itj ; a k , «») - L(u, u j} u k ; a j} a k ) - L(u, u k , Ui\ a k , a*) (3.7) 

taking into account the orientation of the deformation a — > cr', defined as a transition 
between two collections of oriented plaquettes as indicated by the figure. From this 
argument it follows that the independence of the action under such a deformation is 
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locally equivalent to the closure relation fl3.ll) . We consider this invariance an essential 
aspect of the relevant variational principle underlying multidimensionally consistent 
lattice systems. 

The aim of a Lagrangian multiform description over the usual scalar Lagrangian one 
is that it should provide us with not just one variational equation, but in principle an 
arbitrary number of compatible equations. At this stage it is not entirely clear what is 
the optimal formulation for such a principle, in view of the fact that the closure relations 
we have established for specific examples rely on the quadrilateral lattice equations 
themselves, but we will make an attempt in this direction by posing the following. 

Discrete variational principle for integrable lattice systems: The functions u(n) solving 
an integrable multidimensional lattice system on each discrete quadrilateral surface a 
are those for which the action S[u(n); a] of (13. 6p is invariant under local deformations 
of the lattice, as described above, and for which the action attains an extremum under 
infinitesimal local deformations of the dependent variable -u(n). 

The mechanism that we have in mind is as follows. Starting with an action functional 
S[u(n); a] as in (13.61) we impose surface independence of this action. This allows us 
to deform the surface a as we choose, whilst keeping the boundary in place if there 
is a boundary. Thus, we can always render it into a locally flat surface away from 
the boundary, where we can choose any pair of local coordinates ni,nj. In that part 
of the surface we can then apply the usual variational principle with respect to the 
field variables u(n), leading in the usual manner to the Euler-Lagrange equations in 
those lattice directions. If these equations subsequently imply the validity of the closure 
relation for the Lagrangian in terms of which the action is defined, this then ensures that 
those equations are consistent with invariance of the action under deformation of the 
surface which in turn allowed the derivation of those equations in the first place. This 
principle goes farther than just providing a variational derivation of equations of the 
motion from a given Lagrangian, in that in some sense it also imposes conditions on the 
class of admissible Lagrangians to which this principle applies. What is not clear at this 
stage is to what extent admissible Lagrangians can be constructed by application of this 
principle and we do not yet have a general proof that this procedure will automatically 
lead to multidimensionally consistent lattice equations. 

4. Discussion and comparison with the continuous case 

The discrete variational principle formulated in the previous section brings in the 
geometry as a variable of the action functional. This contrasts starkly with the usual 
variational principle, where the Euler-Lagrange equations provide information rather 
on the parametrization of the underlying geometry than on the geometry itself. For 
instance, in the elementary case of a mechanical system with one degree of freedom the 
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action 

rT 



%(*)]= / L(q,q,t)dt (4.1) 



o 

contains hardly any geometry at all, but the relevant Lagrange equation tells us how 
the one-dimensional motion is parametrized in a specific way according to the equations 
of motion. When we have more than one degree of freedom there is obviously room 
for nontrivial phase space geometry, but again the variational equations tell us more 
about how the geometry is parametrized rather than bringing in the geometry as a 
variational variable. Even in classical string theory [13], the geometry of the string 
trajectories (which sweep out a surface in configuration space) plays a role at the level 
of the dependent variables of the string action rather than of the independent variables 
which parametrize the surface. In contrast, our proposal involves the geometry of the 
space of independent variables which is somewhat reminiscent of the de Donder-Weyl 
formalism [IT], although in this approach the connection to integrability is not evident. 
As far as we are aware, all Lagrangian descriptions of (continuous) integrable systems 
so far involve conventional scalar Lagrangians, even where an attempt is made to give 
a multi-Lagrangian description of integrable hierarchies, cf . [35j S3 E] ■ 

At the continuous level, rather than systems of partial difference equations we 
would be looking for systems of partial differential equations. An interesting example 
of a linear system of PDEs which are mutually compatible is given by the following 
nonautonomous set of equations. 

dpAM ~ P 2 j) d K d Pj w = 4 (™A ~ n i d Pj ) -j^-o (njPid Pi - nip)d Vj )w (4.2) 

Pi Pj 

where i,j run over some index set /. Each of these, for fixed labels arise as Euler- 
Lagrange equations from the Lagrange density 



ibjibi \z, Pi y- 



where notably the independent variables Pi,Pj are on an equal footing. Here w = 
w(pi,pj) is the dependent scalar variable and the rij are a pair of parameters of the 
equation, where we associate the parameter with the variable Pi, and the parameter nj 
with the variable pj. It conspires that the system of PDEs (14. 2p . when the labels i,j are 
assumed to run over some index set of cardinality larger than 2, is multidimensionally 
consistent in a similar way as the lattice equations considered in section 2. Furthermore, 
the Lagrangian H4.3f) obeys a closure relation of the following form 

dpXjfc + d Pj L ki + d Pk Lij = (4.4) 

provided one or the other of the following two relations hold 

J n 3 n k p 2 jW Pi n k mp)w P] n t n jP 2 k w Pk \ 

™ m " V (pI - pIM -p)) {pi - p 2 M - pI) (pj -pIM-pDJ [ } 
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or 

(p 2 - p 2 )w PiPj + ( P 2 - p 2 k )w PjPk + (pj-pj)w PkPi = Q 5b 
rtiUj njn k n k ni 

We can infer from the closure relation that the action on the solutions of the system 

S[w; a] — / £j; L jdpi A dpj (4.6) 

is independent of the surface a, relying on Stokes' theorem, and similar to the discrete 
case described in section 3 by locally rectifying the surface, i.e. deforming it locally to 
a plane in terms of selected independent variables Pi,Pj, we can derive from the Euler- 
Lagrange equations in those variables the system of PDEs. It seems somewhat artificial 
in this example to invoke the additional equations (I4.5al) and (14.5b|) . the need for which 



is mainly due to the fact that we are dealing with higher order PDEs in terms of the 
derivatives. We note, however, that the PDEs (14. 2L ( 14. 5 aft and (14.5bl) all hold true on 
a large class of solutions given by the Fourier-type integral of the form 



w 



dk c(k)H(^l) ni (4.7) 



over some suitably chosen curve C in the complex plane and suitably chosen coefficient 
function c(k), where I denotes the index set as above. This example is inspired by the 
canonical form of the plane wave factors, i.e. discrete exponential functions, appearing in 
the solutions of the lattice equations [28], EH] , which explains the use of the notation pi as 
independent variables for historic reasons. The situation described here is the obvious 
continuous analogue of the situation described in section 3 of the relation between 
the closure relation and multidimensional consistency, apart from the fact that we are 
dealing here with a set of linear equations rather than nonlinear ones. 

The full nonlinear case analogous to ( 14.21) appeared first in [30], and it represents 
the full KdV hierarchy as a so-called generating PDE given as follows 

tj tj ( 1 , U Uti UtjtA f i u utj u uu 



n 2 Ul Ul t . i U tt .\ U t . t .U t . t . 

i_ tj i i h_ LlL ] _ I j j 1-1 u " 1 'rj 

"r U UU \ u j. \9 tt2 TT2 + + TT I til i 



(u-ttfuz ui u-tj u ti ) ^ u ti u tj 



n 2 U? U 2 t . l U t . t . 



"\(l,-lj) 2 l? i? t 3 -uu t 

+ 2(t^ tj )* l% {Uti + Utj + 2(tj ~ (4 ' 8) 
n 2 U 1 / 1 1 \ 

+ w^wi [u " + u " + 2(ti - tj)u "" ] + w^-f l » U - ' 

which represents a generalization of the Ernst- Weyl equation of general relativity as was 
shown in [101 . The variables £j are closely related to the Pi of the previous linear example, 
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namely by U — pf . It was argued in [3U| that (14.81) constitutes a multidimensional^ 
consistent system in the same way as the linear equation. The Lagrangian for equation 
( 148]) is 



^=^- t ^+w^Tit +r4 ty (4j) 

This satisfies the closure relation ( 14.4f) provided again that one of two relations hold, 

Uutjtk = 2U U U ftuUtjtkUtjUtkti + UtfltkuUtflJtitj + UtkUtitjUftUtjtk) 

ti i k 

I!? , ™? _j_ ^ C4 10a ) 

2(t k -t i )(t i -t j )Ul 2{t i -t j )(t j -t k )U? J 2(t j -t k )(t k -t i )Ut h ' 

ofl 

{U - t,)U tk U utj + fa - t k )U u U tjtk + (t k - U)U t] U tkU = 0. (4.10b) 

Once again the additional equations f)4.10al) and (|4.10bl) are invoked solely because we 
are dealing with higher order PDEs in terms of the derivatives, which makes it difficult 
to verify by direct computation. Nevertheless all three equations (14. 8p . (I4.10al) and 
(14.10bl) hold on a large class of solutions of soliton type and hence they should certainly 
be compatible between themselves. 



5. Conclusions 



3-point Lagrangians can be identified for all cases in the ABS list, but in this paper 
we have restricted ourselves to those cases for which we have established a closure 
relation. On the basis of this we have formulated a new variational principle in terms of 
Lagrangian multiforms, which we believe captures the multidimensional consistency of 
the underlying integrable systems. Obviously we would want to verify that the closure 
property holds for the remaining cases in the ABS list as well. The cases of Q2 and 
Q3 have a more implicit structure which makes it more difficult to verify the closure 
relation directly. Furthermore the case of the top equation Q4 requires the development 
of new functional identities for the elliptic analogue of the dilogarithm function. We 
intend to deal with those cases in a separate publication. 

An interesting question is whether it is possible to classify integrable discrete 
and continuous systems on the level of the Lagrangians using the closure property. 
Furthermore, we envisage that the formalism proposed in this paper would form 
a paradigm of variational calculus applied to integrable systems, perhaps also in 
connection with associated physical models. In fact, another interesting question is 

Equations (I4.10a|) and (|4.10bj) are manifestations of the fact that 1+1-dimensional equations of KdV 
type can be embedded as dimensional reductions of 2+1-dimcnsional equations of KP type |37j . which 
holds true both for the continuous as well as the discrete case. In fact, the continuous nonautonomous 
equation (|4.10bj) is remarkably similar to the fully discrete Hirota-Miwa equation and we believe that 
it plays the role of a generating PDE for the KP hierarchy. It would be interesting to see how this 
equation fits in with the results of [H?] . 
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what is the quantum analogue of this formalism, e.g. in the context of a path integral 
framework. 
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Appendix A. Dilogarithm functions 

The dilogarithm function is defined by 



Li 2 U) = / ln(1 Z) dz. (A.l) 



Many functional relations involving dilogarithms are given in the book by Lewin[22]. and 
in the review paper of Kirillov[18j, which also covers some of the quantum analogues. 
The pivotal functional relation is the five-term identity 

Li 2 = Li 2 (-^-) + Li 2 (-V-) - U 2 (x) - U 2 (y) 



l — yl — xj \1 — y J \l — x / 

- ln(l - x) ln(l - y), x,y<l. (A. 2) 

For the computations needed for this paper it is more convenient to write (1A.2I) in the 
form 

Li 2 (s) + Li 2 (t) - Li 2 (st) = Li 2 (^T^) + Li2 (~T^i 

+ln (r^) ln (H> s - i>L 

An additional two identities needed are the following, both valid for all real x. 

Li 2 (*) + Li 2 (£) =-i(ln(-x)) 2 -^, (A.4) 

Li 2 (z) + Li 2 (-^_) =-I(ln(l-;r)) 2 . (A.5) 

Equation ( 1A.4I) holds regardless of whether the arguments are positive or negative. 
Equations (IA.2I) and (IA.5[) require additional imaginary terms depending on the sign of 
the arguments; these however cancel out in the course of the closure relation calculations. 
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Appendix B. Proof of the H3 closure relation 

Here we give an outline of the computation needed to show the closure relations (13.11) 
hold for H3. We make use of the dilogarithm identities stated in appendix A. The 
Lagrangian for H3 is 

£aia 2 = &{U, Ux, U 2 , CYi, tt 2 ) 

= -uJ™)+ Li 2 ( ?*L) + Li 2 ( ™) - Li 2 ( ^ 
\— ol\) \—a 2 J \0L\U2) \a 2 u 2 

+ lnQ^ln( M ) + ln(a 2 2 )ln^. (B.l) 

We make a change of variables, similar to those that appear in the 3 leg form of H3. 
This will make the computations simpler and easier to follow. With the abbreviations 

A=m, B=m, C=^ (B.2) 

—OL\ —Ct2 — «3 

the Lagrangian becomes 



Z ai a 2 = -Li 2 (A)+U 2 (B)+U 2 [^)-Li (<X]A 



B ) ~' V a\B 



a\\ 


-A 


1 


-B x 


a\\ 


- B 


1 


-A 2 


a\\ 


-B 


1 


-C2 


a\\ 


-C 


1 


-B 3 


a\l 


-C 


1 


-A 3 



+ In f^j ln(u) + ln(a 2 2 ) In f^j (B.3) 
whilst the equations of motion, written in the variables A,B,C, are as follows: 

(B.4a) 
(B.4b) 

Together with the definitions of A, B, C, these give expressions for A 2 , i?i,etc explicitly 
in terms of A, B, C. To write these in a simple way, define the function Ha,b = 
H(A, B; «i, a 2 ) to be 

_ qj(l - B) - a\{l - A) 
Ha,b (B.5) 

leading to the following 

^3 — —oHha, Bi = —zHab, C 2 = —oHb,Ci 
a{ a 2 a 3 

B C A 

A 2 = —~Habi B3 = -^2~H BC , C\ = —oHca- (B-6) 
af ct 2 cti 

Defining the quantity V as below 

r = a 3 l 

we may now write both the Lagrangians and their shifted versions in terms of A, B and 
C, which leads to 
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Li 2 ( 

V 


S \ . / 

2 Ha,b + Li 2 


ajA\ 
alB ) 


-uJ 

\ 






+ Li 2 


—H B ,c 


+ Li 2 


fa 2 2 B\ 


-Li 2 


—Hb,c] 


+ Li 2 


—He, A 


+ Li 2 


K a\A) 


-Li 2 




) 



+ Li 2 


\aiH AtB J 






-Li 2 




' 2 U B ,c)i 




-Li 2 




-(£ 


>i 







+ In 



..„ v Jln( ffAB ) + lnQ)ln(7/ B , c ) + ln^ i; 

-lng)l„(A)-l„g)ln( B )-l„Q).„( C ) 

— ln(a 2 ) ln(o|) — h^a 2 ,) hi(a|) — ln(cK§) ln(a 2 ) 
+ (ln(a 2 )) 2 + (ln(a 2 )) 2 + (ln(a 2 )) 2 



+ 

where we have we rearranged the terms in a way that suggests which dilogarithm 
identities to use and where. Applying the dilogarithm identity (1A.4|) to the terms 
in the dashed-line boxes, the argument of the dilogarithm functions can be inverted. 
This enables us to use identity (IA.3I) on the terms grouped in the solid-line boxes, using 
the definition of Ha,b to simplify the outcome. We will gather all the logarithm terms 
together at the end. 



+ Li 2 



(A — B)Ha,b 
a\{A-l) 



Lio 



+ Li 2 



[B - C)H B . 



c 



al(B-l) 

(C-A)H c ,a 
af(C-l) 



+ Li 2 
+ Li 2 


(A(B-1)\ 
KB(A-l)) 
(B{C-l)\ 
\C{B-1)) 


r 

!+Li 2 ( 

L 


f C{A-l)\ 
K A{C-l)) 



Lio 



a\{A-l){B-C)H BtC 
al(C-l)(B-A)H A ,B 



Lie 



(C-A)(B-1) 
(B-A)(C-1) 



Li 2 (M-Li/ C - A 



Li; 



(B - A)H A , B 
A(B-C) 



B(A - C) 



- Li 2 



B - A 

a-'b" 

A-C 
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a\{A-l)J \B{A-1)J \a 2 2 {B-l)J \C{B - 1) 
, ,'(4(A-1)\, ( C-A 

, , ' cj(B - 1)(C - A)H C , A \ , 



(B-A)H AB J \B-AJ \B(A-C)J \A-C 
+ ^ + Inf ^ ln(^, fl ) + ln(^f| \n(H B>c ) + In(^) ln(H c>A ) 



f(A- 


1)(B 


-c) 


V(c- 


1)(B 


-A) 


-lnf 


C(A- 


-sy 


B(A - 


-c). 


Hb,c 


))'< 




Hc,a 




■K 1 



6 \ a 2/ V a 3/ 



S) 

KS MA) -KS in(B) -K9 in(c) 



— ln(a^) In (a 2 ;) — ^(a 2 ,) ^M !) — l n ( a 3) hi(ci! 2 ) 

+ (ln(a 2 )) 2 + (ln(a 2 )) 2 + (ln(a 2 )) 2 . (B.8) 

Again, using identity flA.4j) on the terms in the dashed-line boxes, and subsequently 
identity (|A.3[) on the terms grouped in the solid-line boxes, we obtain 



af(C-l) / V af(^-l) 
+ ti ,|M +li /(^* 



al(B-l) J \ al(C-l) 



c4(A-l)J \B{A-1)J \a 2 {B-l)J \C(B-1) 



a 2 3 (C-l)J \A(C-1) 
f al(B-l)(C-A)H C)A \ f (A-l)(B-C) 
\al{C-l){B-A)H A>B ) \{C-l){B-A) 



(B-A)H AB J \B-AJ \B(A-C)J \A-C 

h. i ' c(A . - B i \ m ( {a - m - g ^ + j"j} c - \> i to ( ic - A ^-' 



B(A-C)J \(B-1)(A-C)J \al(A-l)J \(B - A)H AiB 



2V V a\H CtA ) y 2V V ^ca// 2V V C 



2V V C(A-1)JJ 2\ \ (C-A)(B-1)JJ 2\ \ A-B 



+ ln( ln(# AB ) + ln^J ln(tf B ,c) + ln (^ ) HH C ,a) 
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-In 




\n.{a\) 111(02) — l 11 ^) l n ( a 3) — l n ( a 3) m ( a i) 



+ (Hal)) 2 + (\n(al)f + (ln(^)) 2 . 



(B.9) 



Using identity (1A.5j) on the first term of line 1 and the second term of line 2 of (IB. 91) 
leaves the dilogarithm terms which subsequently cancel out. What then remains are 
only the logarithm terms, which also cancel out, leaving L = 0. This concludes the 
proof of the closure relation. 
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